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Abstract: We present a phenomenological study of the associated production of a prompt 
photon and a heavy quark jet (charm or bottom) in Pb-Pb collisions at the LHC. This 
channel allows for estimating the amount of energy loss experienced by the charm and 
bottom quarks propagating in the dense QCD medium produced in those collisions. Cal- 
culations are carried out at next-to-leading order (NLO) accuracy using the BDMPS-Z 
heavy-quark quenching weights. The quenching of the single heavy-quark jet spectrum 
reflects fairly the hierarchy in the heavy quark energy loss assumed in the perturbative 
calculation. On the contrary, the single photon spectrum in heavy-ion collisions is only 
modified at low momenta, for which less heavy-quark jets pass the kinematic cuts. On 
top of single particle spectra, the two-particle final state provides a range of observables 
(photon-jet pair momentum, jet asymmetry, among others) which are studied in detail. 
The comparison of the photon-jet pair momentum, from p-p to Pb-Pb collisions, is sensi- 
tive to the amount of energy lost by the heavy-quarks and could therefore be used in order 
to better understand parton energy loss processes in the heavy quark sector. 

Keywords: heavy ion collisions, parton energy loss, direct photon production, heavy 
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1 Introduction 

Spectacular experimental results on large p T particle production have been reported in Pb- 
Pb collisions at the LHC by the ALICE, ATLAS, and CMS experiments. The quenching of 
light hadrons (inclusive charged hadrons, pions, etc.) has been observed over a very large 
range of transverse momenta [1-3] , soon followed by that of charmed mesons such as D or 
D* [4, 5] . Apart from large-p T hadron production, the reconstruction of jets in heavy-ion 
collisions also allowed for the observation of significant jet asymmetries of jet-jet [6, 7], 
photon-jet [8] and Z-jet [9] correlations, as well as measurements of hadron momentum 
spectra (the so-called fragmentation functions) inside jets [10, 11]. 

Clearly, these results are qualitatively consistent with parton energy loss processes as 
hard quarks and gluons propagate through the dense QCD medium produced in those 
collisions. On the more quantitative side many questions still remain unanswered, among 
which the issue of (radiative) energy loss in the heavy quark sector. Because the heavy 
quark mass acts as a collinear cut-off in the medium-induced gluon radiation, the following 
hierarchy e q > e c > eb of the typical energy loss of light, charm and bottom quarks has 
been proposed by Dokshitzer and Kharzeev [12] and later checked by Armesto, Salgado 
and Wiedemann in the BDMPS-Z energy loss framework [13]. Experimentally, however, 
the slight differences in the quenching of D mesons with respect to that of pions might 
nevertheless be understood as coming from the different color charge of the propagating 
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parton, essentially heavy quarks (Cp = 4/3) and gluons (Ca = 3) for D mesons and pions, 
respectively. As of today, there is no clear and indisputable sign on possible differences 
between the energy loss of light and heavy quarks. On the theoretical side, moreover, the 
above "standard" hierarchy has also been recently questioned by Aurenche and Zakharov 
in Ref. [14] due to formation time effects. It is therefore crucial to identify observables 
at the LHC which could help clarify this issue. In this paper, we argue that measuring 
prompt photon production in association with a heavy quark tagged jet (denoted as 7 + Q 
in the following) in heavy-ion collisions might shed light on the mass dependence of the 
radiative parton energy loss mechanism. 

The production of 7 + Q , first measured in p-p collisions at the Tevatron by the DO 
and CDF experiments [15-18], is a rich and versatile process in various hadronic collisions: 

• In p-p and p-p collisions, first of all, 7 + Q production offers sensitive checks of 
perturbative QCD (pQCD) and might serve as a probe of intrinsic heavy quark 
distributions inside the proton [19]; 

• In p-A collisions (e.g. at RHIC and soon at the LHC), this process can be used 
to constrain the gluon parton distribution function in nuclei (nPDF) [20], which is 
pretty much unconstrained at small values of x. One should underline that knowing 
precisely the nuclear PDFs is a prerequisite in order to obtain reliable predictions in 
heavy-ion collisions; 

• In A-A collisions, finally, the study of 7 + Q provides an ideal tool for investigating 
the energy lost by heavy quarks in the hot QCD medium produced in those collisions. 
Being an electromagnetic probe, the photon produced directly in the hard process 
is expected to traverse the medium unaffected. Its momentum can therefore serve 
as a proxy for the initial momentum of the heavy quark propagating through the 
dense medium and eventually fragmenting into the heavy quark jet. The imbalance 
between the prompt photon and the heavy quark jet momentum going from p-p to 
A-A collisions might thus reflect the amount of energy loss experienced by the heavy 
quark. Furthermore, the comparison between 7 + c and 7 + b production would 
provide access to the mass hierarchy of parton energy loss. 

It is the purpose of this paper to perform an exploratory study of heavy quark energy 
loss from a next-to-leading order (NLO) QCD analysis of 7 + Q production in heavy- 
ion collisions. 1 In addition to the photon and the heavy quark p T spectra, we study the 
distributions in different kinematic variables (qT, Aj, z) which prove useful in order to 
characterize the amount of heavy quark energy loss in the medium. Note that in these 
latter cases, it is mandatory to perform the theoretical calculations at NLO in order to 
obtain reliable results since these variables exhibit a singular behavior at LO accuracy. 

1 A related but different analysis of prompt photon and heavy quark hadrons has been performed pre- 
viously at LO accuracy in the fixed- flavor number scheme [21]. Earlier studies investigating photon corre- 
lations with inclusive hadrons have been performed at NLO, see e.g. [22, 23]. More recently 7 correlations 
with inclusive jets were also investigated [24]. 
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This paper is organized as follows. In Section 2 we describe our theoretical framework 
for the 7 + Q production including the effects of the medium formed during the A-A 
collisions. Our numerical results are presented in Section 3. Finally, in Section 4 we present 
our conclusions and give an outlook on possible future improvements of the analysis. 

2 Theoretical framework 

2.1 7 + Q production in p — p collisions 

Using the standard pQCD collinear factorization framework the cross section for the pro- 
duction of a photon and a heavy quark in p-p collisions is given by 

da = fi ® fj ® dff(ij -»• kQ) <g> D% , (2.1) 

where the sum over all possible partonic subprocesses is performed and fi stands for the 
parton distribution function in the proton. Furthermore, Dl describes the fragmentation 
of the final state parton k into the observed photon and the direct contribution (k = 7) is 
included via D^(z) = 5(1 — z). 

The contributing partonic subprocesses depend on the heavy flavor scheme. The 
present calculation is performed in the variable flavor number scheme (VFNS), which in- 
cludes heavy quark PDFs for factorization scales greater than the threshold given by the 
heavy quark mass, fj,p > rriQ. Furthermore, heavy quark mass terms mg//i 2 have been 
neglected in the perturbative calculation of the partonic cross sections where the factoriza- 
tion/renormalization scale fi is identified with a typical hard scale of the process, that is the 
transverse momentum of the photon or the jet. 2 The range of validity of the calculation is 
thus restricted to momenta \x = O (p T ) larger than a few times the heavy quark mass. Note 
that in heavy-ion collisions the factorization ansatz in Eq. (2.1) is a working assumption. 
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Table 1. List of all direct 2—^3 NLO hard-scattering subprocesses. 

At leading-order accuracy, 0(aa s ), the production of a direct photon with a heavy 
quark jet only arises from the gQ —> 7Q Compton scattering subprocess at the hard- 
scattering level. At next-to-leading order, however, the number of contributing subpro- 
cesses increases to seven (see Table 1). On top of these direct photon subprocesses, the 
contributions coming from the fragmentation of a parton into a photon are also included 
consistently at NLO accuracy. It should however be mentioned that isolation requirements 

2 Note that in the case of a gluon splitting into a collinear QQ pair, the divergency is avoided by imposing 
a cut on the invariant mass such that nigg > 2m,Q [19]. 
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- used experimentally in order to minimize background coming from hadron decays - 
greatly decrease these fragmentation contributions. For further details on the theoretical 
calculations, the reader may refer to [19, 25]. 

2.2 7 + Q production in heavy-ion collisions 

We discuss in this section the calculation of 7 + Q production in heavy-ion collisions, in 
which the final-state partons propagating through the dense QCD medium are expected to 
lose energy through medium-induced gluon radiation, or radiative energy loss processes. 3 
Let us first denote by Pq c and the 4-momenta of the heavy quark and the photon, 
respectively, in p-p collisions. We express these momenta in terms of their transverse 
momentum p™ c and rapidity y vac : 

Pq c = ^(cosh^ ac ,e rQ ,sinh^ ac ) , (2.2) 
= ^(coshy™, e T7 , sinhy™) . (2.3) 

Because of energy loss processes, in heavy-ion collisions the energy of the heavy quark is 
shifted with respect to that in p-p collisions, Eq = E~q c — e, where e is the amount of 
energy lost by the heavy quark while traversing the medium. 4 In the present exploratory 
study, we assume that the photon is not affected by the medium, E-y = E™ c , and moreover 
that the heavy quark does not change its direction while propagating through the medium, 
HQ = 2/g ac . Therefore, the heavy quark momentum in the medium is given by 

PQ = PTQ{coshy Q , e TQ , sinhy Q ) = [p^ - e/ cosh y v Q ac ] (cosh y v Q ac , e TQ , sinhy£ ac ) , (2.4) 

i.e., 

Ptq = Ptq -e T ,e T = e/ cosh y v Q ac . (2.5) 

As discussed in the previous section, beyond the leading order a second parton (labeled 
"2") is produced in real 2 — > 3 subprocesses (see Table 1). The emission of this extra 
parton occurs within the short-distance time scale O [Q ) < 1 fin (where Q 3> Aqcd 
is the scale of the hard process), therefore well before the medium is produced. As a 
consequence, parton 2 also experiences medium-induced energy loss; its energy is thus 
shifted, E2 = E^ — e' , leading to the following expression for its momentum in medium 
similar to (2.4), 

P2 = W2 ~ e7coshyr](coshyr c ,eT2,sinhyr) , (2-6) 
PT2 = p v T a 2 C -4,4 = e7 cosh y 2 vac , V2 = VT , 

after losing the energy e' in the medium. 

3 We do not consider collisional energy loss processes in this study, as they are expected to be less 
important (in comparison to radiative losses) at large momenta. 

4 Strictly speaking, e should be understood as the energy radiated outside the jet cone. We therefore 
implicitly assume in the present study that most of the medium-induced gluon radiation occurs at large 
angles. 
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The momenta of the two final-state partons produced in 2 — > 3 processes are shifted 
independently. In doing so we neglect the possible interference effects in the medium- 
induced gluon radiation off two partons (this question has been addressed, for example, 
in [26]). Also, subsequent (soft and collinear) radiation can nevertheless happen on longer 
time scales and possibly after the radiating parton has escaped the medium (see [27]); such 
radiation is not taken into account in the present fixed-order calculation scheme. Note 
that the energy loss of the extra parton affects the production of 7 + Q events as long 
as it is recombined with the heavy quark Q to form the heavy-flavor jet through a jet 
reconstruction algorithm 5 , and not otherwise. 

In order to compute the medium-modified 7 + Q production cross section in heavy-ion 
collisions, we use the following Monte Carlo procedure: 

1. We obtain an event for the production of a 7 + Q in vacuum at a given heavy quark 



four momentum Pq c and photon four momentum, p. 



vac 
7 ' 



2. For each partonic event (i.e. with one or two partons in the final state), the en- 
ergy loss e (respectively, e') of the heavy quark (respectively, the extra parton) is 
sampled according to a probability distribution or quenching weight, Vi(e), using an 
acceptance- rejection algorithm. 

3. According to Eqs. (2.4) and (2.6), medium-modified four- momenta pq, p 7 and P2 
are constructed. The additional constraints, e < Eq^ c and e' < E^, are moreover 
imposed. 

4. Using the modified four- vectors we evaluate observables {ptq-,Pt^)i as well as the 
correlation variables discussed in Sec. 2.4 in which we bin the events, providing the 
differential cross section for these observables. 

2.3 Quenching weights 

The quenching weights used in the present calculation are obtained using the Poisson 
approximation proposed in [28] , 

»*> - 1 h [n [ ^ A -t\] x * (« - i>) «p {- f £} , m 

n=0 Li=l U U \ i=l / ^ U J 

where j is the flavor of the propagating parton (j = q, g, Q) and d//dw the medium-induced 
gluon spectrum, which in the present work is determined from the perturbative BDMPS-Z 
framework [29]. 

The quenching weights Vj(e) are a scaling function of e/o; c , where ui c = 1/2 q L 2 is the 
typical scale for the energy loss process; the transport coefficient q measures the scattering 
power of the medium (momentum broadening per unit length) and L is the medium length. 
The medium-induced gluon spectrum (and therefore the quenching weights) depends on 
the dimensionless parameter R, R = uj c L, which arises from the kinematic constraint 



5 For the jet reconstruction we combine the heavy quark and parton 2 in a jet if they are within a cone 
of radius R specified in Table 2. 
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restricting the transverse momenta of the radiated gluons. The limit R — >■ oo of the 
spectrum corresponds to the "thick medium" limit of the BDMPS formulation [30] in 
which the opacity, or the number of rescatterings n = L/X — > oo. 6 

In our numerical analysis we use the quenching weights provided in [31] for massless 
quarks and gluons and in [32] for heavy quarks. In the latter case, the weights also depend 
on the heavy quark mass through the dimensionless quantity m/E, where E is the heavy 
quark energy in the medium rest frame. In this exploratory study, we will stick to the 
BDMPS framework, i.e. using the quenching weights computed in the multiple scattering 
scenario (as opposed to the n = 1 opacity case) and taking the "thick" limit, R — > oo. 7 

Due to the heavy quark mass, the medium- induced radiation is suppressed at large 
gluon energies, see [13]. As a consequence, the typical energy loss in this calculation 
is substantially reduced when m/E is large. This calculation confirms in particular the 
hierarchy, e q > e c > e&, conjectured in [12] for heavy quark energy losses in QCD media. 8 As 
mentioned in the introduction, this hierarchy has recently been questioned by Aurenche and 
Zakharov in [14] on the basis of formation time arguments. It is not the goal of the present 
paper to answer this question, nor is it to perform a comprehensive phenomenological study 
which would encompass various assumptions regarding heavy quark energy loss and a full 
treatment of the geometry and the dynamics of the produced medium. Rather we explore 
how the properties of the medium-induced gluon radiation off massive quarks translates 
into a variety of 7 + Q observables through an NLO calculation — within one well defined 
scheme for heavy quark energy loss. 

In the next section, we present the different distributions which we consider in this 
study that we believe would best reflect the dynamics of massive quark energy loss processes 
in heavy-ion collisions. 

2.4 Observables 

The most inclusive observable which should be sensitive to the energy loss is the total 
cross section in dependence of the cut on the ptq (or pr^)- In addition it is useful to 
study the inclusive transverse momentum distributions of the heavy quark and the photon, 
respectively. 

On top of the photon and the heavy quark jet single p T spectra, the two-particle 
final state further offers a wide range of observables which might give a better access to the 
energy loss of the propagating heavy quark. One such variable is the transverse momentum 
difference between the photon and the heavy quark jet, q T , defined as follows: 



Note that in this limit, the probability for no energy loss vanishes. 

r For practical purposes, in order to allow for reasonable CPU time, we use R = 10 5 . This seems to be a 
good approximation since the probability for no energy loss, po, is already below the few percent level for 
the m/E values used in the calculation. 

8 At small values of the heavy quark mass, m/E ~ 10 -2 and intermediate values R ~ 10 3 -10 4 , the energy 
loss of massive quarks is actually larger than in the massive case, questioning the hierarchy e q > e c > et. 
It is however not clear in [13] whether this observation has a well defined physical origin or is an artifact of 
the calculation. 
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At leading order and assuming that the photon is produced directly in the hard process 
(i.e. not by fragmentation), the photon and the jet momenta balance in the transverse 
plane. Therefore this variable reduces at LO to q)j? c = pj£ — = in the vacuum (p-p 
collisions) and g™ ed = pxj ~Ptq = P^y ~~ (.Ptq ~ € t) = e r in the medium (A-A collisions). 
Therefore, a non-zero cross section at a given g™ ed 7^ can be interpreted (at LO, direct) to 
be proportional to the number of events where the heavy quark jet has suffered an energy 
loss ex = g™ ed in the medium. It is unfortunately not as simple because the "vacuum" q T 
distribution in p-p collisions is non-zero at q T > as soon as fragmentation photons or 
NLO contributions are taken into account (the former being however suppressed from the 
use of isolation criteria) . 

Assuming the photon and the jet to be in different hemispheres (see also the cuts in 
Table 2) the momentum difference in A-A collisions becomes q™ ed = p^ — (p^q — £t) = 
q™ c +€T- In other words the q T spectrum in A-A collisions is shifted by -Kt as compared to 
p-p collisions, where ey should be understood as a typical amount of (transverse) energy 
loss. The shift should therefore be more pronounced for light jets than for charm and 
bottom quark jets, because of the hierarchy used in the calculation. 

Similar to the photon-jet pair momentum we shall investigate as well the imbalance 

Aj, 

Aj = P^-PTQ , (2.9) 

PT-y + PTQ 

which has been measured for dijet [6, 7, 33] and photon-jet correlations [8] by ATLAS and 
CMS. Finally the distribution in the momentum imbalance, 

z lQ -2 (2.10) 

PTy 



is also considered. In the leading order kinematics (and assuming the photon is produced 
directly in the hard process), these two variables would reduce respectively to Af ed = 

„vac , \ „„J ~mcd -i , /^vac 



ex/(2p^ — ex) and z^q = 1 — er/Pr^ in heavy- ion collisions. As compared to the q T 
variable, these are sensitive to the ratio of the energy lost over the parton energy (or, 
fractional energy loss), er/p^, instead of the absolute magnitude of the energy loss, ex- 
Generically, for all the above observables quenching factors -Raa ar e defined as 

1 dq(A + A^ 7 + Q) 

-KAA = -ir T — -T-, . — , (2.11) 

iVcoii dcr(p + p -> 7 + Q) 

where N co \\ stands for the number of binary nucleon-nucleon collisions (N co \\ = A 2 in 
minimum bias A-A collisions). In particular, we consider the suppression factors for the 
inclusive pT-y spectrum, E? AA (pT-y), and for the inclusive ptq spectrum, RjIa(Ptq)- 

Some comments are in order: (i) With this definition the ratios Raa an d Rj^a would 
in general be different from unity even in the absence of cold and hot nuclear matter effects 
due to differences between parton distributions inside protons and neutrons as is the case 
for inclusive photon production (see e.g. [34]). However, in our case the dominant initial 
state is g + Q and these parton distributions should be (to a very precise level) the same 
inside protons and neutrons, (ii) To construct this ratio experimentally, the p-p cross 
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section which is measured at a different center-of-mass energy has to be extrapolated to 
the center-of-mass energy of the A-A run. For this reason it is extremely valuable that the 
photon transverse momentum spectrum which is insensitive to hot medium effects (up to a 
sensitivity on the cut on the heavy quark jet energy) can be used to calibrate the effects of 
the energy extrapolation and of cold nuclear matter, (iii) An advantage of considering the 
double ratios R^2 = ^aa/^aa 15 * na ^ uncertainties due to the choice of scales and nuclear 
PDFs should largely cancel. While the scale uncertainties already cancel to a good degree 
in the simple ratios there remains a residual uncertainty due to the lack of knowledge of 
the nuclear gluon distribution, see e.g., [20]. Hence, it will be crucial to use data from the 
upcoming p-Pb run at the LHC at the beginning of 2013 to better constrain the nuclear 
gluon distribution. 9 Therefore, for the time being we refrain from using nuclear PDFs 
[35-38] in this study where we focus on the effects of the heavy quark energy loss solely. 

3 Results 

3.1 Ingredients 

The present calculations have been carried out using the CTEQ6.6M PDFs inside a proton 
[39] (along with the corresponding strong coupling constant a s ' 5 (Mz) = 0.118 at next-to- 
leading order) and the photon fragmentation functions of Bourhis, Fontannaz and Guillet 
[40]. The renormalization, factorization and fragmentation scales have been set to [Ir = 

= fif = px-y and we have used m c = 1.3 GeV and m,b = 4.2 GeV for the charm and 
bottom quark masses. All cross sections (for both p-p and Pb-Pb collisions) have been 
calculated at a center-of-mass energy ^/s NN = 5.5 TeV. 

For our exploratory numerical studies, we apply the kinematic cuts provided in Table 2. 
Most notably, the cuts on the minimal pt^ and ptq are asymmetric in order to avoid 
configurations with ~ ptq where the NLO cross section is known to become infrared 
sensitive [20]. For simplicity, we consider events at mid-rapidity in small bins around y 7 = 
and HQ = 0. Moreover we impose a cut on the azimuthal angle between the photon and 
heavy quark jet, A^ 7 q > 37r/4 in order to ensure that the photon and the heavy quark are 
produced in different hemispheres. Once data are available with sufficient statistics these 
results could be updated with the appropriate experimental cuts. 
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Rapidity 
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Jet radius 


Photon 
Heavy quark jet 


p™™ = 20 GeV 
p4% = 12 GeV 


|y 7 | < 0.2 
IvqI < 0-2 


R = 0.4, e < 0.1.E7 7 


R = 0A 



Table 2. Phase space cuts used for the theoretical predictions in our study. 



Alternatively, constraints on the nPDFs could possibly be obtained through inclusive photon and weak 
boson production in A-A collisions. 
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Table 3. Total integrated cross section and number of events per year for 7 + Q production in 
Pb-Pb collisions at the LHC. 



3.2 Total cross sections and event rates 

In this section we compute the NLO cross sections for p+p — > 7+c+X and p+p — > 7+6+X 
taking into account the cuts on the phase space listed in Table 2. The corresponding cross 
sections for Pb-Pb collisions are obtained by scaling the ones for p-p collisions with A 2 
where A = 207 is the atomic number of lead. In Table 3, we present results for the cross 
section in vacuum and for two media with parameters R = 10 5 , u c = 50 GeV and R = 10 5 , 
co c = 100 GeV, respectively. The expected number of events per year are given in the 
third column assuming a yearly luminosity of £pKp b = 0.5 nb -1 for Pb-Pb collisions at the 
LHC [41]. As can be seen, the event numbers are sufficiently large to measure the cross 
section for both, 7 + c and 7 + 6 production where the 7 + 6 cross section is roughly a factor 
of 7 smaller. The smaller cross sections in medium are due to the fact that fewer events 
pass the cut ptq > Ptq when energy loss is present. With sufficiently large statistics, it 
could also be interesting to analyze the cross section as a function of p™g • The 7 + c cross 
section clearly depends on uj c . The suppression factor of the total cross section is given by 
Raa = 0.87 (.Raa = 0.74) for uj c = 50 GeV (co c = 100 GeV). Conversely, the dependence 
of the 7 + 6 cross section on uj c is much more modest with Raa = 0.95 (Raa = 0.93) 
for uj c = 50 GeV (co c = 100 GeV). It is noteworthy that the suppression effect could be 
further enhanced by optimizing the values for the cuts on pt^ and ptq- In particular, the 
quenching could be enhanced (without reducing the cross section too much) by imposing 
an upper integration bound pt^ < p™j X ^ 50 GeV as can be inferred from Fig. 1(a). 

The numbers in Table 3 can be used to obtain rough estimates of cross sections and 
event rates when other phase space cuts are employed or detection efficiencies for the 
photon and the heavy quark are taken into account: (i) Assuming an approximately flat 
rapidity dependence, the cross section will roughly scale with the size of the rapidity bins 
of the photon and the heavy quark jet. We checked that using |y 7 | < 2 and \uq\ < 2 leads 
to a cross section roughly 60 times larger than the one quoted in Table 3; (ii) Assuming a 
power-law behavior for the px-y and ptq dependence (see Sec. 3.3.2) the cross sections will 
depend strongly on the values oip™ 1 ™ and Ptq- Results for p™ 1 ™ and p^q different from 
the ones in Table 2 can be estimated by multiplying the cross sections/event numbers with 
a factor (12 GeV/p^) n x (20 GeV/p^ 7 n ) n with n ~ 4. 

While a more precise estimate of the expected event numbers clearly depends on the 
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Figure 1. Next-to-leading order differential cross section a) da/dpx-y as a function of px 1 and b) 
da/dpxQ as a function oipxQ in pb/GeV. Shown are results for the vacuum (black solid lines) and 
for different media with R = 10 5 and w c = 50 (red dashed lines) and 100 GeV (blue dash-dotted 
lines) . 



details of the experimental acceptances and efficiencies our results indicate that the number 
of 7+c events produced in a year of Pb-Pb collisions at the LHC will be substantial and that 
also a measurement of the 7 + 6 cross section will be feasible before the projected upgrade 
of the luminosity during the second long shutdown of the LHC where an instantaneous 
luminosity £p b p b = 6 x 10 27 cm -2 s _1 is expected [42] corresponding to a yearly (t = 10 6 s) 
luminosity of £p e ^p b = 6 nb _1 . 

3.3 Single-inclusive transverse momentum spectra 
3.3.1 Differential cross sections 

Before turning to the nuclear production ratios, we first present the predictions for the 
single p T spectra in p-p and Pb-Pb collisions. In Fig. 1 the NLO differential cross sections 
da/dpT-y (left) and da/dpTQ (right) for 7 + c production are shown for the vacuum and 
for different media (co c = 50 and 100 GeV, with R = 10 5 ). As discussed in Sec. 2, in our 
framework the photon is not affected by the medium. However, as can be seen in Fig. 1(a), 
at small pt~/, the pr-y spectrum in medium is reduced below the vacuum spectrum. The 
reason is that fewer events pass the cut on the heavy quark jet transverse momentum 
{ptq > 12 GeV), if heavy quark energy loss is present. This is a physical effect which 
scales with the heavy quark energy loss at small PTy'- at lo c = 50 GeV (oo c = 100 GeV) we 
observe a 15% (30%) reduction. Conversely R P ^2 * s unity at values away from the cut, say 
for pt~/ > 70 GeV, since for these momenta most of the heavy quark jets pass the cuts as 

PTQ >PTQ- 

The ptq spectra are much more affected by the medium over the entire px range, as 
is clearly visible in Fig. 1(b). In this case, the quenching is not related to the presence of 
the heavy quark threshold required in the calculation but rather due to the shift of the jet 



- 10 - 




80 100 120 140 160 180 

p Tl [GcV] 



100 150 

Vtq [GcV] 



(a) 



(b) 



Figure 2. Same as in Fig. 1 for the 7 + 6 case. 



momentum, see Eq. (2.4), in heavy-ion collisions. At pxq < 50 GeV the heavy quark mass 
mildens the energy loss effects, leading to a small difference between the vacuum cross 
section (black solid line) and the in medium cross sections (red dashed line, blue dash- 
dotted line) . This will be further discussed in the next section where quenching factors are 
presented. 

In Fig. 2 we present the same observables, but for 7 + b production. As is visible the 
same trends are present as in Fig. 1, which are however clearly reduced in size compared 
to the 7 + c case. 

3.3.2 Quenching factors Raa 

We consider here the quenching factor R^a^Ptq) °f * ne heavy-quark jet. The goal is 
twofold. Firstly, as mentioned in Section 2.3, the aim is to show the sensitivity of the 
observables involving a photon and heavy quark jet in the final state on parton energy 
losses entering via different (massless and massive) quenching weights. Secondly, we want 
to explore the effects of the NLO QCD corrections on the quenching factor R® a (ptq)- 

As an illustration of the first point, the quenching of the charm jet p T spectrum 
in 7 + c production has been computed as a function of pt c hi Fig. 3(a) under various 
hypothesis regarding the energy loss of charm quarks, i.e. assuming that the charm quark 
suffers the same energy loss as that of a light quark (blue dash-dotted line) or that of a 
heavy quark assuming mq = m c (red dashed line) and mq = m& (black solid line) in the 
quenching weights (for completeness the calculation has also been carried out using the 
quenching weight of a propagating gluon, brown dotted line). The calculation is performed 
at leading order direct, i.e. according to the Compton subprocess only. The quenching 
factors computed using these various prescriptions follow closely the hierarchy of parton 
energy loss in the BDMPS-Z quenching weights used for each flavor in the calculation, which 
are shown in Fig. 3(b). In particular, a much stronger (respectively, weaker) charm jet 
quenching would be expected should the charm quark lose energy like a gluon (respectively, 
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bottom quark). Interesting also is the comparison between the predictions assuming light 
quark and charm quark energy loss. At low transverse momenta, the charm quark energy 
loss predicts less suppression as compared to the light quark energy loss, while at large 
momenta, pt c 3> m c, the induced energy loss of a charm and a light quark become identical, 
as can be seen in Fig. 3(a) (by comparing the dashed and dash-dotted lines). 10 



7+c (Leading Order) 

_R = 10 5 , w =50 GcV 



100 



(a) 




(b) 



Figure 3. a) Suppression factors Raa in leading order versus transverse charm momentum, pt c , 
applying energy loss according to the quenching weights shown in the right figure; b) Quenching 
weights for energy loss corresponding to gluons (dotted, brown line), light quarks (dash-dotted, 
blue line), charm quarks (dashed, red line), and bottom quarks (solid, black line) in dependence of 
the energy loss e scaled to io c . 



It is possible to extract information on the parton energy loss from the suppression 
factors .Raa in Fig. 3(a), by assuming a power law dependence of the differential cross 
section da/dpTQ vs. ptq- 

da™ c „ , , 

where the index n increases smoothly from n ~ 4 to n ~ 5.5 from low to high transverse 
momenta. The medium cross section is obtained by shifting the transverse momentum by 
a typical transverse energy loss ey (which could have a ptq dependence) : 

j med Hrr vac 

( PT Q med ) * ? (PTQ = PTQ mCd + Z T ) ■ (3-2) 



dp T Q med ' dpTQ 

As a consequence the suppression factor is approximately given by 

RQ AA(PT Q )=(-^y-) n , (3.3) 

\PTQ + £T / 

leading to 

Zt = Ptq ((Ri A r l/n - l) , Z = ~e T cosh y Q . (3.4) 



Similarly, the solid black line merges with the quenching factor for light quarks at even larger pr c 
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Table 4. Extracted energy loss (column 4) from the nuclear suppression factor Raa (column 2) 
in Fig. 3(a) using Eq. (3.4) for bottom, charm, light quark, and gluon energy loss, at two different 
Ptq values (column 1). For comparison the peak energy loss (column 6) and the mean energy loss 
(column 7) scaled to lu c — 50 GeV are shown as well, see text for details. 



Using Eq. (3.4) and the values for .Raa in Fig. 3(a) at ptq = 30 GeV (taking n = 4) 
and ptq = 90 GeV (n = 5.5) we obtain the energy loss values e presented in column four 11 
of Table 4 (in the fifth column the value is scaled by w c = 50 GeV). We compare this 
to the energy loss corresponding to the peak of the quenching weights for bottom, charm, 
light quark and gluons in column 6. Finally in the seventh column we show the energy loss 
corresponding to the mean of these quenching weights. 

It is interesting to note that the ratio of the energy loss for gluon and light quarks 
(13.5/6.0 or 13.0/6.0) is very close to the ratio of the color factors Ca/Cf = 9/4. As 
can be seen the energy loss for light quarks and gluons is independent of the chosen ptq 
value, whereas the heavy quarks lose less energy at small Ptq due to the mass effects. 
Furthermore at ptq = 90 GeV the charm quark energy loss is already almost identical to 
the light quark one. As can be seen the extracted values for the energy loss (in column five) 
are systematically larger (by a factor varying between 1.4 and 3.0) than the most likely 
( "peak" ) energy loss according to the quenching weights (in column 6) and systematically 
smaller (by a factor varying between 1.2 and 1.6) than the average energy loss. 

Now turning to the NLO results, we present predictions for the quenching factors Raa 
for both, 7 + c (red dashed and solid lines) and 7 + 6 (blue dotted and dash-dotted lines) 
production, in Fig. 4(a). It is reassuring to see that the quenching factors at NLO are not 
very different from the ones obtained at LO (compare, for example, the red dashed lines in 
Figs. 4(a) and 3(a)). The corresponding double ratios -Raa/-^aa f° r w c = 50 (dashed line) 
and 100 GeV (solid line) are shown in Fig. 4(b). Clearly, the double ratio for uj c = 100 GeV 
is much smaller than the one for co c = 50 GeV implying a strong sensitivity to the energy 
loss of the charm and bottom quark jets in the medium. Interestingly, the two double 
ratios have a very similar shape for both values of ui c . 
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Figure 5. Next-to- leading order differential cross section da/dqT for the production of 7 + c for a) 
negative and b) positive values of qr- Shown are results for the vacuum (black solid lines) and for 
different media with R = 10 5 and w c = 50 (red dashed lines) and 100 GeV (blue dash-dotted lines). 
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Figure 6. Feynman diagrams for the subprocesses a) gQ —> gQj and b) gg —> QQj giving the 
main contribution to the qr distributions in Figs. 5 and 8 at positive values of qx- 
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Figure 7. Feynman diagrams for the subprocesses a) gg — > QQ7 and b) — > 5Q7 and — >• 
giving the main contribution to the gy distributions in Figs. 5 and 8 at negative values of qT- 



3.4 Two-particle observables 

We now turn to a discussion of the two-particle observables introduced in Sec. 2.4. We 
begin with the qx distribution defined in Eq. (2.8) for 7 + c production which is shown in 
Fig. 5(a) for negative and in Fig. 5(b) for positive values of qx- We have spared the region 
—20 GeV < qT < 20 GeV where pt 1 — Ptq to avoid problems with IR-safety 12 as discussed 
in Sec. 3.1. In addition to the vacuum cross section (black solid line) results are shown for a 
medium with u c = 50 GeV (red dashed line) respectively ui c = 100 GeV (blue dash-dotted 
line) . As can be seen the medium spectra are right-shifted with respect to the vacuum one 
for both, negative and positive values of qx (see the discussion in Sec. 2.4). Clearly, the shift 
increases with increasing oj c and is smaller for positive qx where it is about 2 GeV for oj c = 
50 GeV and varying between 3.5 and 5 GeV for u c = 100 GeV. For negative qT the shift is 
roughly 5 GeV (11-15 GeV) for lo c = 50 GeV (u c = 100 GeV). The overall differences in 
shape and size of the curves in Figs. 5(a) and 5(b) can be understood as a consequence of 
the different subprocesses which contribute (dominantly) to the cross section for positive 
and negative qT- For positive qx, the dominant contribution comes from the subprocess 
gQ — > gQ'j followed by gg —> QQj where appropriate Feynman diagrams are shown in 
Figs. 6(a) and 6(b), respectively, where the leading photon is balanced by the Qg pair in 
Fig. 6(a) and by the QQ pair in Fig. 6(b). For negative qT, the dominant contribution 
comes from the subprocess gg — > QQj where in this case the relevant Feynman diagram is 
shown in Fig. 7(a) for a leading heavy quark jet. 13 Other non- negligible subprocesses (for 
qT < 0) are gQ — > gQj and qQ — > qQj (see Fig. 7(b) for sample diagrams). 

The larger shift in energy for negative qx (compared to the positive qx case) can 
be understood in the following way: the photon momentum pj-y remains fixed close to 
its minimal value ]?T7 mm = 20 GeV whereas ptq varies (i.e. is larger) in order to satisfy 
q T ~ p Tl mm —ptq. Hence the heavy quark energy loss is less pronounced since mass effects 
are less important at larger ptq values. Conversely, for q? > 0, the situation is inverted 
where now ptq remains fixed close to its minimal value PTQ min = 12 GeV whereas px-y 
varies (i.e. is larger) in order to satisfy qx — pr-y — PTQ mm - In this case, the heavy quark 

11 Note that e ~ It for central rapidities \y\ < 0.2. 

12 Note that meaningful results can still be obtained when integrating/binning over a sufficiently large 
region around qr = as has been done in Fig. 9 for the observables Aj (binning around Aj — 0) and z 
(binning around 2=1). 

13 Note that this amplitude is suppressed in the case of positive qr because we require that the tagged 
heavy quark be the one with the largest pr- 
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Figure 8. Same as in Fig. 5 for 7 + 6 production. 



energy loss will be smaller because the heavy quark mass effects are still relevant. In Fig. 8 
we present the corresponding results for 7 + 6 production. Again, the same trends can be 
observed which are, however, reduced in size. 

Finally, in Fig. 9 we show results for the jet asymmetry Aj (left) and the momentum 
imbalance z (right). The variable Aj defined in Eq. (2.9) is very similar to qr and the 
Aj spectrum shows indeed very similar features as the qx distribution, in particular the 
in- medium curves are shifted to the right, the shift is larger for negative Aj, and the 
distribution is slightly asymmetric around Aj = 0. The z-distribution peaks at z = 1 
which corresponds to a configuration where the heavy quark and the photon are back-to- 
back (as in LO). With our cut on the angle between the heavy quark jet and the photon 
momentum (6 > 37r/4) we find 0.7 x < z = —^2- cos 9 < Therefore, roughly, 

the region z < 1 (z > 1) corresponds to ptq < pr-y (ptq > Pt^)- In this case, the in- 
medium spectra are left-shifted and the bigger shifts are visible in the region z > 1 which 
corresponds, as in the previous figures, to the kinematic configuration ptq > pr-y- 

4 Conclusions 

An important aspect of jet quenching phenomenology is to understand the mass dependence 
of parton energy loss. While the quenching of single inclusive hadron production might give 
a clue on this topic, the possibilities offered by the LHC experiments allow for investigating 
more exclusive processes which have the advantage that a range of interesting observables 
can be constructed from the momenta of the final state particles. 

In this paper, we performed an exploratory study of the associated production of 
a photon and a heavy quark jet in heavy-ion collisions, at NLO accuracy. On general 
grounds this is a very promising process in which the energy loss of the heavy quark jet 
can be calibrated using the prompt photon momentum since the latter is assumed to be 
unaffected by the hot and dense QCD medium. Despite the smaller cross sections as 
compared to the single inclusive heavy-quark jet production, the expected rate at the LHC 
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Figure 9. Differential cross section with respect to a) the jet asymmetry Aj and b) the momentum 
imbalance z int the 7 + c production channel. 



is sufficiently large to ensure reasonable statistical uncertainties for both 7 + c and 7 + 6 
production. In addition to the inclusive pr spectra of the photon and the heavy quark 
jet, we have performed NLO calculations for a range of observables (qT, Aj, z) in order to 
identify which of these distributions provide useful information on the amount of energy 
loss experienced by the propagating heavy quark. 

The inclusive ptq and pT-y distributions are very promising observables. Numerical 
results for da/dpT-y and da/dpTQ calculated at NLO QCD were presented. We find that 
the ptq spectra are affected by the medium over the entire pr range, whereas the px-f 
distribution only gets quenched at small px-y due to the cut on the jet transverse momentum. 
Comparing these two spectra should thus be particularly instructive in order to disentangle 
the effects due to the heavy-quark energy loss from the ones due to the nuclear modification 
of the parton distribution functions and other cold nuclear matter effects. Information on 
the parton energy loss could be obtained either by studying the total cross section as a 
function of the cut on the jet transverse momentum or by analyzing the ptq distribution. 
For the latter we have presented in Fig. 4(a) the quenching factors Raa for both, 7 + c and 
7 + 6 production at NLO QCD. Very interesting is also the double ratio R%x/^AA shown 
in Fig. 4(b) of which the shape (not the normalization) as function of ptq = PTc = PTb 
turns out to be almost independent of the choice of the parameter uj c . 

The distributions in the two-particle kinematic variables have also been investigated in 
detail, paying attention in particular to the photon-jet pair momentum q T . As expected, 
these distributions are shifted towards larger values in heavy-ion collisions. The comparison 
of q T spectra in p-p and Pb-Pb collisions for 7 + c and 7 + 6 production should thus allow 
for a "direct" access to the amount of energy lost by charm and bottom quarks, respec- 
tively. Note that such distributions are singular at leading order accuracy and therefore 
NLO predictions prove mandatory in order to compute such spectra at large (positive and 
negative) q T values. 



-17- 



For completeness the distributions in the photon-jet asymmetry, Aj, and momentum 
imbalance, z, have also been determined in the 7 + c channel in p-p and Pb-Pb collisions. 
Similar patterns to the ones observed in the q T distributions are reported. The comparison 
between the distributions in various kinematic variables should thus help to determine at 
a quantitative level the amount of heavy-quark energy loss from the future high-precision 
measurements through detailed phenomenological studies. 

Let us mention some future improvements to be carried out beyond the present study. 
The systematic comparison of 7 + c/b jet production with 7+ inclusive jet production (for 
which first measurements have been reported recently [8]) should be particularly inter- 
esting and useful. In addition, as soon as the nuclear gluon parton distribution is better 
constrained from data of the forthcoming p-A run at the LHC, the predictions should 
be updated using nuclear PDFs including systematic studies of the nPDF and scale un- 
certainties. Finally, it would also be important to compare the present calculations with 
calculations based on other energy loss frameworks and different assumptions regarding 
the modeling of the heavy-quark energy loss in the medium. 

Acknowledgment 

We are grateful to thank S. Peigne for useful discussions on parton energy loss. This 
work was supported by the CNRS through a PICS research grant and the Theorie-LHC- 
France initiative. This work is funded by "Agence Nationale de la Recherche" under grant 
ANR-PARTONPROP. 



- 18 - 



References 

[1] ALICE Collaboration, K. Aamodt et. al., Suppression of Charged Particle Production at 
Large Transverse Momentum in Central Pb-Pb Collisions at y?s NN = 2.76 TeV, Phys. Lett. 
B696 (2011) 30, [1012.1004]. 

[2] ATLAS Collaboration, Measurement of the centrality dependence of charged particle spectra 
and rep in lead-lead collisions at snn = 2. 76 tev with the atlas detector at the Ihc, Tech. Rep. 
ATLAS-CONF-2011-079, CERN, Geneva, Jun, 2011. 

[3] CMS Collaboration, S. Chatrchyan et. al.. Study of high-px charged particle suppression in 
Pb-Pb compared to p-p collisions at ^fs^n =2.76 TeV, Eur. Phys. J. C72 (2012) 1945, 
[1202.2554]. 

[4] ALICE Collaboration Collaboration, Z. Conesa del Valle, D mesons suppression in PbPb 
collisions at ^snn = 2.76 TeV measured by ALICE, 1210.2163. 

[5] ATLAS Collaboration, Measurement of the centrality dependence of open heavy flavour 
production in lead-lead collisions at sqrts = 2. 76 tev with the atlas detector, Tech. Rep. 
ATLAS-CONF-2012-050, CERN, Geneva, May, 2012. 

[6] ATLAS Collaboration, G. Aad, Observation of a Centrality- Dependent Dijet Asymmetry in 
Lead-Lead Collisions at -^/snn = 2.76 TeV with the ATLAS Detector at the LHC, Phys. Rev. 
Lett. 105 (2010) 252303, [1011.6182]. 

[7] CMS Collaboration, S. Chatrchyan et. al., Observation and studies of jet quenching in PbPb 
collisions at nucleon-nucleon center- of -mass energy = 2.76 TeV, Phys. Rev. C84 (2011) 
024906, [1102.1957]. 

[8] CMS Collaboration, S. Chatrchyan et. al., Studies of jet quenching using 

isolated-photon+jet correlations in PbPb and pp collisions at -^/saT/v = 2.76 TeV, 1205.0206. 

[9] ATLAS Collaboration, Measurement of momentum imbalance in z —¥ ££ + jet events in 
lead-lead collisions at i/snn— 6 2- 76 tev with the atlas detector, Tech. Rep. 
ATLAS-CONF-2012-119, CERN, Geneva, Aug, 2012. 

[10] CMS Collaboration Collaboration, S. Chatrchyan et. al., Measurement of jet 

fragmentation into charged particles in pp and PbPb collisions at sqrt(s[NN])= 2.76 TeV, 
JEEP 1210 (2012) 087, [1205.5872]. 

[11] ATLAS Collaboration, Measurement of inclusive jet charged particle fragmentation 

functions in pb+pb collisions at sqrt(s(nn)) = 2.76 tev with the atlas detector, Tech. Rep. 
ATLAS-CONF-2012-115, CERN, Geneva, Aug, 2012. 

[12] Y. L. Dokshitzcr and D. Kharzeev, Heavy quark colorimetry of QCD matter, Phys. Lett. 
B519 (2001) 199-206, [hep-ph/0106202]. 

[13] N. Armesto, C. A. Salgado, and U. A. Wiedemann, Medium- induced gluon radiation off 
massive quarks fills the dead cone, Phys. Rev. D 69 (2004) 114003. 

[14] P. Aurenche and B. Zakharov, Anomalous mass dependence of radiative quark energy loss in 
a finite-size quark- gluon plasma, JETP Lett. 90 (2009) 237-243, [0907.1918]. 

[15] DO Collaboration, V. M. Abazov et. al, Measurement of 7 + b + X and 7 + c + X 
production cross sections in pp collisions at y/s=1.96 TeV, Phys. Rev. Lett. 102 (2009) 
192002, [0901.0739]. 



- 19 - 



[16] CDF Collaboration, T. Aaltonen et. al., A Study of the associated production of photons and 
b-quarkjets in pp collisions at y/s = 1.96 TeV, Phys.Rev. D81 (2010) 052006, [0912.3453]. 

[17] DO Collaboration, V. M. Abazov et. al., Measurement of the photon-\-b-jet production 

differential cross section in pp collisions at ^fs — 1.96 TeV, Phys.Lett. B714 (2012) 32-39, 
[1203.5865]. 

[18] DO Collaboration, V. M. Abazov et. al, Measurement of differential photon+ c-jet cross 
sections and the ratio of photon+ c and photon+ b cross sections in proton- antiproton 
collisions at *Js = 1.96 TeV, 1210.5033. 

[19] T. Stavreva and J. Owens, Direct Photon Production in Association With A Heavy Quark At 
Hadron Colliders, Phys.Rev. D79 (2009) 054017, [0901.3791]. 

[20] T. Stavreva, I. Schienbein, F. Arleo, K. Kovarik, F. Olness, et. al., Probing gluon and 

heavy-quark nuclear PDFs with 7 + Q production in pA collisions, JHEP 1101 (2011) 152, 
[1012.1178]. 

[21] Z.-B. Kang and I. Vitev, Photon-tagged heavy meson production in high energy nuclear 
collisions, Phys.Rev. D84 (2011) 014034, [1106.1493]. 

[22] F. Arleo, P. Aurenche, Z. Belghobsi, and J. -P. Guillet, Photon tagged correlations in heavy 
ion collisions, JHEP 11 (2004) 009, [hep-ph/0410088]. 

[23] H. Zhang, J. F. Owens, E. Wang, and X.-N. Wang, Tomography of high-energy nuclear 

collisions with photon- hadron correlations, Phys. Rev. Lett. 103 (2009) 032302, [0902.4000]. 

[24] W. Dai, I. Vitev, and B.-W. Zhang, Momentum imbalance of isolated photon-tagged jet 
production at RHIC and LHC, 1207.5177. 

[25] T. Stavreva, Direct photon production in association with a heavy quark. PhD thesis, Florida 
State University, 2009. URN: etd-04102009-155714. 

[26] Y. Mehtar-Tani, C. A. Salgado, and K. Tywoniuk, Anti-angular ordering of gluon radiation 
in QCD media, Phys. Rev. Lett. 106 (2011) 122002, [1009.2965]. 

[27] J. Casalderrey-Solana, J. G. Milhano, and P. Q. Arias, Out of Medium Fragmentation from 
Long-Lived Jet Showers, Phys.Lett. B710 (2012) 175-181, [1111.0310]. 

[28] R. Baier, Y. L. Dokshitzer, A. H. Mueller, and D. Schiff, Quenching of hadron spectra in 
media, JHEP 0109 (2001) 033, [hep-ph/0 106347]. 

[29] U. A. Wiedemann, Gluon radiation off hard quarks in a nuclear environment: Opacity 
expansion, Nucl.Phys. B588 (2000) 303-344, [hep-ph/0005129]. 

[30] R. Baier, Y. L. Dokshitzer, A. H. Mueller, S. Peigne, and D. Schiff, Radiative energy loss of 
high-energy quarks and gluons in a finite volume quark - gluon plasma, Nucl.Phys. B483 
(1997) 291-320, [hep-ph/9607355]. 

[31] C. A. Salgado and U. A. Wiedemann, Calculating quenching weights, Phys.Rev. D68 (2003) 
014008, [hep-ph/0302184]. 

[32] N. Armesto, A. Dainese, C. A. Salgado, and U. A. Wiedemann, Testing the color charge and 
mass dependence of parton energy loss with heavy-to-light ratios at RHIC and CERN LHC, 
Phys.Rev. D71 (2005) 054027, [hep-ph/0501225]. 

[33] CMS Collaboration, S. Chatrchyan et. al., Jet momentum dependence of jet quenching in 
PbPb collisions at ^fs^ = 2.76 TeV, Phys.Lett. B712 (2012) 176-197, [1202.5022]. 



- 20 - 



[34] F. Arleo, Hard pion and prompt photon at RHIC, from single to double inclusive production, 
JEEP 09 (2006) 015, [hep-ph/0601075]. 

[35] K. J. Eskola, H. Paukkunen, and C. A. Salgado, EPS09 - A New Generation of NLO and LO 
Nuclear Parton Distribution Functions, JEEP 04 (2009) 065, [0902.4154]. 

[36] M. Hirai, S. Kumano, and T. H. Nagai, Determination of nuclear parton distribution 

functions and their uncertainties at next-to-leading order, Phys. Rev. C76 (2007) 065207, 
[0709.3038]. 

[37] I. Schienbein, J. Yu, K. Kovarik, C. Keppel, J. Morfin, et. ai, PDF Nuclear Corrections for 
Charged and Neutral Current Processes, Phys.Rev. D80 (2009) 094004, [0907.2357]. 

[38] D. de Florian, R. Sassot, P. Zurita, and M. Stratmann, Global Analysis of Nuclear Parton 
Distributions, Phys.Rev. D85 (2012) 074028, [1112.6324]. 

[39] P. M. Nadolsky, H.-L. Lai, Q.-H. Cao, J. Huston, J. Pumplin, et. ai, Implications of CTEQ 
global analysis for collider observables, Phys.Rev. D78 (2008) 013004, [0802.0007]. 

[40] L. Bourhis, M. Fontannaz, and J. Guillet, Quarks and gluon fragmentation functions into 
photons, Eur. Phys. J. C2 (1998) 529-537, [hep-ph/9704447] . 

[41] F. Arleo et. al, Photon physics in heavy ion collisions at the LEC, hep-ph/0311131. 

[42] ALICE Collaboration, B. Abelev et al, "Upgrade of the ALICE Experiment. Letter Of 
Intent." http://cdsweb.cern.ch/record/1475243/files/LHCC-I-022.pdf, 2012. 



- 21 - 



